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Building  integrated  photovoltaic  (BIPV)  systems  are  a  relevant  application  of  photovoltaics.  In  countries 
belonging  to  the  International  Energy  Agency  countries,  24%  of  total  installed  PV  power  corresponds  to  BIPV 
systems.  Electricity  losses  caused  by  shadows  over  the  PV  generator  have  a  significant  impact  on  the 
performance  of  BIPV  systems,  being  the  major  source  of  electricity  losses.  This  paper  presents  a  methodology 
to  estimate  electricity  produced  by  BIPV  systems  which  incorporates  a  model  for  shading  losses.  The 
proposed  methodology  has  been  validated  on  a  one  year  study  with  real  data  from  two  similar  PV  systems 
placed  on  the  south  facade  of  a  building  belonging  to  the  Technical  University  of  Madrid.  This  study  has 
covered  all  weather  conditions:  clear,  partially  overcast  and  fully  overcast  sky.  Results  of  this  study  are  shown 
at  different  time  scales,  resulting  that  the  errors  committed  by  the  best  performing  model  are  below  1%  and 
3%  in  annual  and  daily  electricity  estimation.  The  use  of  models  which  account  for  the  reduced  performance 
at  low  irradiance  levels  also  improves  the  estimation  of  generated  electricity. 
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1.  Introduction 

Over  the  past  two  decades  the  number  of  grid  connected  PV 
systems  has  experienced  an  important  growth.  According  to  the 
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International  Energy  Agency  IEA  [1],  between  1992  and  2010  the 
average  growing  rate  of  installed  PV  power  was  39%  in  its  member 
countries.  This  growth  has  been  due  to  the  inherent  benefits  of 
renewable  energies  and  supporting  polices  by  public  entities, 
mostly  based  on  feed-in-tariffs.  The  total  cumulative  power 
installed  in  IEA  countries  at  the  end  of  2010  was  35  GW.  Of  this 
capacity,  nearly  24%  was  integrated  in  or  associated  to  buildings. 

Despite  progressive  reductions  of  feed-in-tariffs  in  most  coun¬ 
tries,  it  is  expected  for  BIPV  systems  to  keep  their  growth  rate. 
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With  rising  prices  of  retail  electricity  and  decreasing  costs  of  PV, 
grid  parity  with  commercial  electricity  will  soon  become  a  reality, 
especially  in  urban  residential  areas.  This  fact,  together  with  new 
policies  that  enhance  self-consumption  (currently  under  discussion 
in  Spain  [2])  points  that  building  integrated  photovoltaic  (BIPV) 
systems  will  become  more  relevant  in  the  future.  Also  to  be 
mentioned  is  the  new  European  Energy  Performance  of  Buildings 
Directive  [3],  where  PV  can  definitely  play  a  new  role  in  the  so- 
called  “nearly  zero-energy  buildings”  which  are  described  in  the 
directive  as  buildings  with  a  “very  high  energy  performance.  [...]  The 
nearly  zero  or  very  low  amount  of  energy  required  should  be  covered  to 
a  very  significant  extent  by  energy  from  renewable  sources,  including 
energy  from  renewable  sources  produced  on-site  or  nearby". 

The  increase  in  installed  PV  power  demands  for  simple, 
accurate  models  able  to  predict  the  electricity  expected  from  a 
system.  These  models  should  be  confronted  with  actual  data  from 
existing  installations  to  validate  them.  However,  in  the  literature 
there  is  a  lack  of  specific  studies  about  performance  of  PV  systems 
integrated  in  buildings  [4-6],  As  a  consequence  there  is  also  a  lack 
of  models  which  describe  system  behaviour  with  reasonable 
accuracy. 

Location  determines  the  performance  of  a  PV  system  in  two 
ways:  climate  and  obstacles  surrounding  the  system.  Available 
irradiance  and  ambient  temperature  are  set  by  climate  while 
obstacles  may  reduce  the  effective  irradiance  over  the  PV  gen¬ 
erator.  Orientation  and  slope  of  the  PV  generator  surface  also  affect 
the  effective  irradiance  received.  However,  while  energy  loss 
due  to  small  deviations  on  orientation  and  slope  from  optimal 
position  have  low  sensitivities  [7],  the  impact  of  shading  on 
system  performance  is  relevant. 

Shadows  over  BIPV  systems  are  casted  by  nearby  objects  such 
as  trees,  utility  poles,  surrounding  buildings  or  even  the  building 
itself  on  which  the  system  is  integrated.  Regardless  their  nature, 
these  obstacles  are  characterized  by  irregular  shape  and  closeness 
to  the  PV  generator.  The  irregular  shape  makes  its  modelling  and 
incorporation  into  the  obstacle  profile  a  complex  and  tedious  task. 
The  closeness  of  obstacles  to  the  PV  generator  means  that  an 
important  fraction  of  the  sky  dome  is  occupied  by  these  obstacles 
and  shading  has  a  high  impact  on  BIPV  systems.  Therefore, 
a  practical  estimation  method  of  BIPV  systems  performance 
should  include  an  accurate  methodology  for  evaluating  shading 
losses.  This  paper  presents  a  methodology  to  estimate  the  elec¬ 
tricity  produced  by  a  BIPV  system  based  on  the  combination  of  a 
novel  model  for  the  estimation  of  shading  losses  with  well-known 
validated  models  for  estimation  of  electricity  production.  The 
performance  of  this  methodology  has  been  assessed  with  a  case 
study:  The  hourly  electricity  produced  over  one  year  by  two  BIPV 
systems  integrated  in  the  same  building,  one  heavily  shaded  while 
the  other  is  relatively  free  of  shadows,  has  been  estimated  using 
the  methodology  and  compared  to  the  actual,  measured  values. 

The  structure  of  the  remainder  of  the  paper  is  as  follows.  Section 
2  describes  the  procedure  to  estimate  system  output,  paying  special 
attention  to  the  model  for  shading  losses.  Section  3  describes  a  BIPV 
system  used  to  validate  this  methodology.  Results  and  comments  on 
this  validation  are  given  in  Section  4.  Finally,  conclusions  on  the 
results  and  methodology  usage  are  given  in  Section  5. 

I.J.  Background 

Any  method  for  the  estimation  of  shading  losses  is  based  on  the 
combination  of  two  steps  that  can  be  considered  as  two  indepen¬ 
dent  processes: 

-  First,  the  irradiance  over  the  surface  of  interest  -  the  PV 

generator  -  is  obtained.  This  step  needs  as  input  two  sets  of 

data:  a  series  of  annual  irradiation,  typically  hourly  values,  and 


a  description  of  the  obstacles  surrounding  the  PV  generator. 
The  comparison  of  obstacle  profile  and  sunpath  determines  at 
which  hours  the  generator  is  obscured  by  shadows. 

-  On  the  second  step  the  electricity  that  can  be  extracted  from 
the  PV  generator  is  estimated.  This  is  expressed  in  terms  of  a 
dimensionless  factor  called  shading  loss  factor:  LSh  which  can 
be  added  to  the  various  losses  also  existing  in  a  PV  system. 

On  the  past  decades  several  methods  developed  to  address 
shading  losses  have  been  published  on  the  literature.  Some  of 
these  methods  include  the  two  steps  described  above  or  are 
limited  to  only  one  of  the  steps. 

Alonso-Garda  [8,9],  Di  Pizza  [10],  Karatepe  [11],  Kawamura  [12] 
and  Rauschenbauch  [13]  provide  theoretical  models  to  simulate 
the  behaviour  of  PV  modules  and  generators  in  conditions  of  non- 
uniform  radiation.  However,  they  do  not  provide  methods 
to  obtain  the  distribution  of  shadows  over  the  surface  under 
study.  The  majority  of  these  models,  the  only  exception  being 
Kawamura's  one,  take  into  account  the  effect  of  bypass  diodes. 
Also,  three  authors  (Rauschenbauch,  Kawamura  and  Alonso- 
Garda)  include  the  effect  of  the  inverse  /-V  characteristic  of  the 
shaded  cell  to  estimate  the  1-V  overall  curve  of  the  module  or 
generated  affected  by  shadows. 

Other  authors,  on  the  contrary,  provide  models  to  estimate  the 
effective  irradiation,  that  is,  after  incorporating  shading  effects  but 
these  models  do  not  estimate  the  electricity  losses  in  which  the  PV 
system  incurs  due  the  effect  fn  inhomogeneous  irradiation.  Drif 
[14]  proposes  a  model  based  on  geometrical  description  of 
obstacles  in  which  considers  each  component  of  solar  radiation, 
beam  and  diffuse,  separately.  Drif's  model  has  been  validated 
experimentally  on  a  PV  system  located  at  Jaen  University  (Spain). 
Niewenda  [15]  introduces  a  software  tool,  SOMBRERO  that  allows 
to  obtain  the  irradiation  over  any  surface  in  the  presence  of 
obstacles.  This  tool  is  not  limited  exclusively  to  PV  systems  and 
can  be  combined  with  other  tools  ( SUNCODE ,  TRNSYS)  to  evaluate 
the  impact  of  shadows  on  passive  applications  of  solar  energy. 

Finally,  Woyte  [16]  and  Oozeki  [17]  analyze  shading  losses  from 
the  perspective  of  electricity  production  of  the  overall  PV  system, 
focusing  their  analysis  on  electricity  losses  rather  than  irradiation 
losses.  Woyte  studies  the  impact  of  shadows  on  a  PV  generator 
according  to  its  electrical  configuration:  generator  formed  by 
several  strings  connected  in  parallel  to  one  central  inverter,  each 
string  connected  to  an  inverter  or  AC  modules.  For  large  obstacles 
whose  shadow  covers  a  significant  section  of  the  generator  Woyte 
does  not  find  differences  in  performance  between  the  three 
configurations.  Oozeki  introduces  a  model  to  estimate  electric 
losses  caused  by  shadows,  as  well  as  other  losses,  based  on  the 
monitoring  of  the  system  during  one  year.  Oozeki's  model  great 
advantage  is  that  it  does  not  require  a  characterization  of  the 
obstacles.  On  the  other  hand,  its  main  drawback  is  that  it  requires 
a  time  span  of  one  year  and  it  cannot  be  used  to  estimated  shading 
losses  during  the  design  phase  of  a  PV  system. 

2.  Estimation  of  shading  losses 

Fig.  1  shows  an  outline  of  the  procedure  proposed  to  estimate 
PV  system  output.  Shading  losses  are  the  most  relevant  factor  on 
BIPV  systems  performance  and,  therefore,  the  model  to  estimate 
them  is  the  core  of  the  methodology  shown  in  the  figure.  The 
model  for  the  estimation  of  shading  losses  has  been  developed 
according  to  the  following  objectives: 

-  Accuracy:  obstacle  profile,  despite  its  complexity,  must  be 
characterized  accurately  by  the  methodology;  this  will  result 
in  better  estimation  of  associated  losses. 
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Fig.  1.  Flowchart  of  methodology  for  system  output  estimation. 


_  Simplicity:  the  methodology  must  be  flexible  and  simple 
enough  to  describe  obstacles,  however  complex  they  are,  using 
few  and  easy-to-understand  information. 

The  models  in  Fig.  1  are  described  in  detail  below. 

2.1.  Solar  radiation  model 

The  three  components  of  solar  irradiation  over  an  arbitrarily 
oriented  and  tilted  surface  can  be  estimated  if  the  global  irradia¬ 
tion  over  the  horizontal  plane  is  known.  For  this  purpose  scheme 
proposed  by  Liu  and  Jordan  [18]  is  used:  beam  and  diffuse 
components  of  horizontal  irradiation  are  estimated  from  global 
irradiation  through  the  existing  correlation  between  clearness 
index  and  diffuse  fraction.  The  use  of  a  local  correlation,  when 
available,  is  preferred  over  the  use  of  a  global  correlation  [7,19], 
Beam  and  albedo  irradiation  components  over  an  inclined 
surface,  B^a.p)  and  Rt{a,p)  respectively,  during  a  given  time  period 
t  are  estimated  in  a  straightforward  manner  as  shown  in  Eqs. 
(1)  and  (2),  where  Bt( 0)  is  beam  component  of  horizontal  irradia¬ 
tion,  Gt(0)  is  global  component  of  horizontal  irradiation,  6S  is 
incidence  angle,  8ZS  is  solar  zenith  angle,  a  is  orientation  angle 
(0  for  south  orientation  and  positive  to  the  west),  p  is  tilt  angle 
(0°  for  horizontal  surface)  and  p  is  surface  albedo: 

Bt(a,p)  =  Bt( 0)  ■  max(0,  cos  8S)/  cos  9ZS  (1) 

Rt(a,P)  =  Gt(0)  •  (1  —  COS  P)  ■  p/2  (2) 

Perez  model  is  used  to  estimate  diffuse  irradiation  [19],  This 
model  has  been  selected  over  similar  models  due  to  its  better 
performance  [7,20-22], 

The  model  described  so  far  needs  global  irradiation  over  a 
horizontal  surface  to  estimate  the  components  of  irradiation  over  a 
tilted  surface  which  is,  in  turn,  necessary  to  estimate  shading 
losses.  However,  in  B1PV  systems,  if  irradiation  is  measured  it  is 
always  available  on  the  same  orientation  and  tilt  as  the  PV 
generator.  Therefore,  a  procedure  to  estimate  horizontal  irradia¬ 
tion  from  in-plane  irradiation  is  necessary  since  there  are  no 
models  in  the  literature  for  the  correlation  between  diffuse  and 


global  components  of  irradiation  over  tilted  surfaces.  A  novel 
approach  to  obtain  values  of  horizontal  irradiation  based  on 
selecting  days  from  a  typical  meteorological  year  (TMY)  from  the 
PV  system  location  is  proposed.  This  approach  has  the  advantage 
that  it  relies  on  data  sets  that  characterize  statistically  the  local 
climate  conditions  and  are  increasingly  available  [23,24],  Every 
day  of  measured  tilted  irradiation  is  replaced  by  a  day  from  the 
TMY.  It  is  not  possible  to  replace  an  actual  day  with  the  same  day  - 
same  date  -  from  the  TMY  because  the  diffuse  fraction  of  both  days 
will  not  be  the  same.  Instead,  every  day  is  replaced  by  the  day 
from  the  TMY  which  is  most  similar  to  it.  Similarity  between  days 
is  estimated  by  the  function  in  Eq.  (3 ),/s.  The  first  term  inside  the 
brackets,  ptmy.data .  is  the  correlation  coefficient  between  both 
days;  this  term  assures  that  the  days  from  the  TMY  will  have  a 
similar  energy  distribution  than  the  days  from  the  recorded  set. 
While  correlation  coefficient  provides  days  with  similar  power 
profiles,  it  does  not  guarantee  that  these  days  have  the  same 
irradiation.  For  this  purpose  a  comparison  between  daily  energies 
is  included  as  second  term;  Gd  TMY(0)  is  daily  irradiation  of  the  day 
from  the  TMY,  Gddata{a,p)  is  daily  irradiation  of  the  actual  day  and  II 
II  represents  absolute  value  function.  The  search  inside  the  TMY  is 
restricted  to  a  window  with  a  size  of  2w+l  days,  centred  on  the 
day  to  be  replaced  and  extended  w  days  ahead  and  behind  of  it. 
From  within  this  interval  the  day  that  maximizes  fs  is  selected: 

fs  =  IPtmy.data  +  U  —  W(GijM\(G)  —  Gd  iata(a,  p))  /  Gd  data(a,  p)\\\) /2  (3) 

2.2.  Shading  model 

The  proposed  methodology  for  shading  losses  estimation 
is  based  on  the  following  concepts: 

-  Irradiation  components  are  considered  separately.  Global  irra¬ 
diation  is  separated  into  its  three  components:  beam,  diffuse 
and  ground-reflected. 

-  Obstacles  and  PV  generator  are  modelled  as  three  dimensional 
objects.  On  B1PV  systems  the  size  and  distance  of  the  obstacle 
to  the  generator  are  frequently  similar  to  the  dimensions  of  the 
generator  itself.  Therefore,  it  cannot  be  assumed  that  the 
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Fig.  2.  Grid  defined  over  the  PV  generator  surface  for  shading  losses  estimation  and 
translation  procedure  of  obstacle  coordinates  between  points  inside  the  PV 
generator. 


generator  is  constrained  to  a  single  point;  instead,  the  gen¬ 
erator  is  considered  as  a  2  dimensional  surface.  In  addition,  the 
PV  generator  will  not  be  evenly  lighted  when  shaded:  some 
parts  will  be  covered  by  shadows  while  others  will  be  illumi¬ 
nated  directly  by  the  Sun.  In  contrast,  the  much  simpler 
assumption  that  the  generator  is  dimensionless  only  allows 
for  two  possible  states:  fully  illuminated  or  fully  covered  by 
shadows. 

-  Shading  losses  are  estimated  using  geometrical  considerations. 
Although  the  effect  of  non-uniform  lighting  conditions  on  the 
electrical  behaviour  of  a  PV  system  is  not  linear  with  the 
surface  fraction  covered  by  shadows  [8-13]  the  model  devel¬ 
oped  only  takes  into  account  the  geometrical  configuration  of 
the  PV  generator,  not  its  electrical  configuration.  The  use  of  an 
electrical  model  to  account  for  shading  losses  in  a  PV  system 
requires,  in  addition  to  estimate  parameters  associated  to  the 
model,  that  physical  and  electrical  configuration  of  both  PV 
modules  (number  of  cells,  series-parallel  configuration,  size 
and  number  cells  protected  by  a  single  bypass  diode)  and 
generator  (series-parallel  configuration,  placement  of  modules 
along  their  long  or  short  side).  This  situation  can  be  compli¬ 
cated  even  further  if  PV  modules  are  placed  following  an 
irregular  pattern,  as  it  is  the  case  on  some  BIPV  systems  for 
constructive  or  aesthetic  purposes.  On  the  contrary,  the  geo¬ 
metric  approach  proposed  provides  a  less  complex  shading 
model:  the  only  information  needed  concerning  the  PV  gen¬ 
erator  is  its  physical  dimensions,  orientation  angle  (a)  and  tilt 
angle  (/?). 


As  shown  in  Fig.  1,  shading  losses  are  estimated  following 
a  two-step  process.  Both  steps  are  described  below. 


2.2.1.  First  step:  irradiation  over  PV  generator  surface 

First  step  goal  is  to  estimate  the  incident  irradiation  on  any 
point  of  the  PV  generator  surface.  This  step  starts  by  defining  a 
grid  of  points  over  the  PV  generator.  Obstacles  and  generator  sizes 
are  assumed  to  be  of  similar  magnitude,  meaning  that  every  point 
of  the  generator  “sees"  a  different  obstacle  profile.  The  points  in 
the  grid  are  selected  so  that  they  are  coincident  with  the  corners  of 
the  modules  in  the  generator  (as  indicated  by  the  grey  dots  in 
Fig.  2).  For  every  point  in  the  grid  an  associated  obstacle  profile  is 
calculated.  An  obstacle  profile  is  a  two  dimensional  function  that 
assigns  to  every  point  in  the  sky  a  binary  value:  0  if  that  point  is 
occupied  by  an  obstacle,  1  otherwise.  The  range  of  this  function 
comprises  values  of  azimuth  from  - 180°  to  180°  and  elevation 
from  0°  to  90°.  The  obstacle  profile  is  defined  as  discrete  function 
so  it  can  be  easily  computed  by  software.  The  resolution  is  set  at 
0.5°  for  azimuth  and  0.25°  for  elevation.  The  obstacle  profile 
function  is  formally  defined  in  Eq.  (4),  where  i//  is  azimuth  and  y 


is  elevation. 


OP(i//,r) 


0,  3  obstacle 

1  otherwise 


\ji  =  [  —  180, 180],  y  =  [0. 90) 
(4) 


Fig.  2  shows  the  procedure  to  obtain  the  obstacle  profile  for 
every  point  in  the  grid.  Initially,  the  coordinates  of  all  relevant 
obstacle  points  relative  to  one  point  of  the  PV  generator  surface,  Gj 
in  Fig.  2,  are  known.  To  obtain  the  obstacle  profile  seen  from  point 
G2  all  points  in  the  obstacle  are  translated  from  the  coordinate 
system  centred  in  Gi  to  a  new  coordinate  system  whose  origin  is 
G2.  The  translation  vector  is  the  coordinates  of  G ,  with  respect  to 
G2;  vector  G2Gi.  For  example,  the  coordinates  of  point  O j,  which 
belongs  to  the  obstacle,  with  respect  to  G2  (r2)  are  obtained  by 
adding  the  translation  vector  to  the  coordinates  of  Ch  relative  to 
Ci: 

r2  =  f*i +  G2Gj  (5) 


Once  all  points  of  the  obstacle  have  been  translated,  their 
rectangular  coordinates  are  converted  into  spherical  coordinates: 
elevation  ( y )  and  azimuth  (i/r).  The  repetition  of  this  translation 
process  for  all  the  points  in  the  grid  provides  a  set  of  obstacle 
profiles.  This  set  of  obstacle  profiles  contains  all  the  necessary 
information  to  know  the  distribution  of  shadows  over  the 
PV  generator  surface  at  any  moment. 


2.2.2.  Second  step:  estimation  of  shading  loss  factor 

The  second  step  estimates  losses  due  to  shading  on  a  PV 
generator  from  the  set  of  obstacle  profiles  generated  in  the  previous 
step.  It  starts  by  evaluating  individually  the  effect  of  obstacles  on 
eveiy  component  of  solar  irradiation:  beam  B,  albedo  R,  circumsolar 
diffuse  Dcs,  isotropic  diffuse  Dls  and  horizon  band  diffuse  DH.  Each 
component  has  a  different  behaviour  regarding  obstacles,  by  evalu¬ 
ating  them  separately  five  partial  shading  loss  factors  are  obtained. 

-  Beam  and  circumsolar  diffuse  irradiation  share  in  common 
that  they  come  from  the  Sun  position  in  the  sky.  Therefore, 
their  respective  shading  loss  factors,  LSh,b  and  LSh,cs.  are 
estimated  in  the  same  way.  The  grid  superimposed  to  the 
PV  generator  divides  its  surface  into  rectangular  regions.  The 
behaviour  of  each  region  is  governed  by  the  four  points  at  its 
vertexes:  if  one  region  has  two  or  more  of  its  vertexes  in  the 
shaded  part  of  the  PV  generator  this  region  does  not  con¬ 
tribute  to  the  electricity  production  of  the  generator.  The 
partial  shading  loss  factor  for  beam  and  circumsolar  diffuse 
is  then  defined  as 

(number  of  shaded  regions} 

sh.b  —  sh.cs  —  jtota]  number  of  regions)  * 


The  reason  why  a  shaded  region  does  not  affect  the  behaviour 
of  non-shaded  regions  is  because  it  is  assumed  that  bypass 
diodes  effectively  isolate  the  effect  that  shaded  modules  have 
on  non-shaded  modules  of  the  PV  generator.  By  identifying 
eveiy  module  with  a  region  it  is  assumed  that  only  one  bypass 
diode  is  installed  per  module  that  limit  the  electrical  power 
losses  to  the  ratio  between  the  area  protected  by  one  single 
diode  and  total  PV  module  area  [25,26],  This  is  a  conservative 
approach  since  modules  usually  have  several  bypass  diodes 
installed.  Flowever,  with  this  assumption  it  is  not  necessaiy  to 
know  how  many  bypass  diodes  a  module  has,  thus  reducing 
the  complexity  of  shading  effects  calculation  [27], 

-  Isotropic  diffuse  irradiation  comes  from  all  points  of  sly 
hemisphere.  For  eveiy  point  in  the  grid  its  associated  loss 
factor  lSHjs  is  the  ratio  between  the  solid  angle  of  the  obstacles 
and  the  total  solid  angle  of  the  sky  hemisphere  seen  by  the  PV 
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generator.  The  subscript  in  1Sh,is  denotes  that  this  loss  factor  is 
defined  for  a  single  point  not  for  the  whole  PV  generator.  It 
must  be  noted  that  the  amount  of  sky  seen  by  a  PV  generator 
depends  on  its  tilt  angle.  It  is  considered  for  convenience  that, 
in  the  estimation  of  the  solid  angle  occupied  by  obstacles,  the 
radial  distance  is  equal  to  1.  Therefore,  it  is  omitted  in  Eq.  (7). 
The  shading  loss  factor  for  the  whole  PV  generator,  LSFIS,  is  the 
average  of  all  partial  factors  obtained  from  the  PV  generator 
points  defined  by  the  grid: 

I  (solid  angle  of  obstacles)  Xi//Xf{OP(i//,  y)  ■  sin  y/  ■  Ayr  ■  Ay) 

SFJS  ~  (solid  angle  of  sky  hemisphere)  —  2 In-fi) 

(7) 

Lsfjs  =  mean(/SFB)  (8) 


-  Horizon  band  diffuse  and  albedo,  although  originated  from 
different  physical  phenomena,  have  the  same  behaviour:  they 
both  come  from  the  horizon  line.  The  partial  loss  factors  (Ish.h: 
Ish,r)  for  every  point  in  the  grid  are  defined  as  the  fraction  of 
the  horizon  seen  by  the  PV  generator.  A  PV  generator  always 
sees  half  of  the  horizon  line,  an  arc  of  n  radians  centred  on 
generator  orientation,  not  depending  on  its  inclination  and 
orientation.  As  before,  the  loss  factors,  LSf,h  and  LSFr,  are  the 
averages  of  the  corresponding  partial  loss  factors: 


,  ,  (angle  of  obstacles) 

sf,h  -  SFjt  -  (ang|e  0f  horizon) 


=  £lf( OP(v,r)  A¥)  with 
n 


l,y  =  0 

(9) 


f-sF,H  =  mean(/SFH);LSFR  =  mean(/SFR) 


(10) 


Finally,  the  overall  shading  loss  factor,  LSF ,  is  obtained  from  the 
weighted  combination  of  the  loss  factors  estimated  for  every 
solar  irradiation  component: 

f-SF  =  (fit  •  isF.B  +  Dtcs  •  Lsf.cd  +  Dus  ■  LSFjs+Dch 

■  LSF.H+Rt  ■  LSf,r)/G  (11) 


2.3.  Power  models 

The  electricity  produced  by  a  PV  system  can  be  determined  by 
estimating  first  the  maximum  DC  electricity  that  can  be  delivered 
by  the  PV  generator  and  later  subtracting  losses  related  to  different 
physical  phenomena: 

-  Shading  losses  Lsh-  these  losses  are  described  in  Section  3.2. 

-  Optical  losses  Lco:  losses  caused  by  the  reflection  of  incident 
light  when  incidence  angle  is  different  than  perpendicular.  The 
model  proposed  by  Martin  has  been  used  [28], 

-  Electrical  losses  LCF.  electrical  losses  include  losses  due  to 
voltage  drops  in  wires  on  the  DC  side  and  PV  modules 
mismatch.  In  this  study  electrical  losses  are  included  in  the 
models  for  maximum  DC  electricity  estimation  and  they  are  not 
considered  separately. 

-  System  losses  Ls:  these  losses  account  for  the  instantaneous 
DC-AC  conversion  efficiency  of  the  inverter.  The  widely  used 
model  proposed  by  Schmidt  has  been  selected  [29]  because  its 
theoretical-practical  approach,  which  allows  to  estimate  mod¬ 
el's  parameters  from  actual  measurements  of  inverter  opera¬ 
tion.  This  model  is  presented  in  the  following  equation: 

V  =  Po/Pi  =  Po/(Po  +  losses)  =  po/(po  +  k0  +  ki  po+k2p20 )  (12) 

where  ;;  is  inverter  efficiency,  pa  is  normalized  inverter  output 
power,  Pi  is  normalized  inverter  input  power  and  k,  are 


coefficients  associated  to  conversion  losses:  ka  represents 
inverter  self  consumption,  k,  represents  inverter  losses  whose 
dependence  on  power  is  linear  (diode  consumption,  switching 
devices)  and  k2  represents  losses  whose  dependence  on  power 
is  quadratic  (wires,  resistors,  coils,  etc.). 

-  Thermal  losses:  since  the  effect  of  temperature  on  system  yield 
is  included  in  the  models  for  maximum  DC  electricity  estima¬ 
tion  thermal  losses  are  not  accounted  independently. 

AC  electricity  produced  by  a  PV  system  is  estimated  by  the 
following  equation,  which  incorporates  the  losses  previously 
defined: 

Eac  =  Edc,max(G,Tc)  ■  (1  — Lsh )  •  (1  —Lco)  ■  (1  —Lce)  ■  (1  —is)  (13) 

Maximum  DC  electricity  produced  by  a  PV  generator,  Eoc.max , 
depends  only  on  operation  conditions:  irradiance  and  cell  tem¬ 
perature.  In  this  study  two  different  models  for  the  estimation  of 
maximum  DC  electricity  have  been  tested.  The  first  model  was 
proposed  by  Osterwald  [30],  This  model  combines  simplicity  and 
accuracy  and  performs  better  than  more  complex  models  for 
crystalline  silicon  modules  [31],  Osterwald  model  is  presented  in 
the  following  equation: 

Edc,max  =  fine  ■  Gt  •  [1  -yp  ■  ( Tc-Tc*)]/G *  (14) 

where  Gt  is  irradiation,  G*  is  irradiance  at  standard  test 
conditions  (1000  W/m2),  Tc  is  cell  temperature,  Tc*  is  cell  tem¬ 
perature  at  standard  test  conditions  (25  °C),  PnC  is  generator  peak 
power  at  standard  test  conditions  and  yP  is  temperature  coefficient 
of  maximum  DC  power. 

Osterwald  model  has  been  used  in  combination  with  the 
nominal  operation  cell  temperature  model  NOCT  [32],  presented 
in  the  following  equation: 

Tc  =  TA  +  (NOCT-20)Gt/800  (15) 

where  TA  is  ambient  temperature  and  NOCT  is  nominal  opera¬ 
tion  cell  temperature.  The  model  based  on  NOCT  parameter  has 
been  selected  because  its  performance  is  similar  to  more  complex 
models  for  estimating  cell  temperature  [33],  In  addition,  this 
parameter  is  available  from  the  PV  modules  datasheets. 

During  this  study,  Osterwald  model  was  found  to  be  inaccurate 
for  low  levels  of  irradiance.  In  consequence,  a  model  that  defines 
PV  module  efficiency  as  a  function  of  irradiance  was  included  in 
the  study.  The  expression  of  Eq.  (16),  first  proposed  by  Beyer  in 
[34],  has  been  selected,  where  e,  are  fitting  parameters  that  must 
be  evaluated  for  every  PV  system.  The  effect  of  cell  temperature  is 
the  same  proposed  by  Osterwald. 

0  =  >l(Gt,  Tc)  =  [ei  +e2  ■  Gt+e3  •  log ( Gt)]  •  [1  —  y  ■  (Tc  —  Tc*)]  (16) 


3.  Experimental  validation 

3.1.  PV  system  description 

The  methodology  proposed  in  Section  2  has  been  validated 
using  a  BIPV  system  located  at  the  Technical  University  of  Madrid 
(“Lopez  Araujo"  building,  ETSI  Telecomunicacion).  The  modules 
are  placed  above  the  windows  of  the  building  south  facade  acting 
as  sunshades;  in  winter,  when  the  sun  is  low  on  the  sky,  the 
sunlight  shines  through  the  windows  helping  to  warm  up  the 
building;  in  summer,  with  the  sun  higher  in  the  sky,  the  PV 
generator  blocks  direct  sunlight  into  the  building  reducing  thermal 
load.  Fig.  3  shows  this  system. 

The  system  is  formed  by  282  PV  modules,  oriented  8“  East  and 
tilted  32°,  distributed  in  13  strings  of  serial  connected  modules. 
Nine  of  these  strings  are  connected  in  parallel  to  a  central  inverter 
and  the  four  remaining  strings  are  connected  to  string  inverters. 
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Fig.  3.  PV  system  located  on  "Edificio  Lopez  Araujo"  building  at  ETS1.  Telecomu- 
nicacion  Madrid.  (For  interpretation  of  the  references  to  colour  in  this  figure 
caption,  the  reader  is  referred  to  the  web  version  of  this  paper.) 


The  system  also  includes  a  meteorological  station  with  the 
following  sensors:  two  reference  PV  sensors  of  the  same  technol¬ 
ogy  as  the  PV  generator  (crystalline  silicon)  for  the  measurement 
of  irradiance  [35]'1  and  cell  temperature  [36]  and  a  PT-100  resistor 
for  the  measurement  of  ambient  temperature.  Both  PV  sensors  are 
of  same  technology  as  the  modules,  with  same  orientation  and  tilt 
angles  as  the  modules  and  located  on  the  roof  top  on  a  place  free 
of  shadows. 

For  this  study  information  from  two  of  the  individual  strings, 
SI  and  S2,  composed  of  26  series  connected  PV  modules  each 
(outlined  in  red  in  Fig.  3),  as  well  as  data  from  the  meteorological 
station,  has  been  collected  and  analyzed.  The  reason  for  selecting 
these  strings  is  that  they  are  representative  of  BIPV  systems  where 
the  PV  generators  act  as  sun  shading  elements,  an  architectural 
function  that  often  implies  partial  shading  in  the  PV  modules.  In 
fact,  PV  modules  of  string  S2  are  shaded  during  afternoon  hours  by 
the  upper  row  SI.  Thus,  obstacle  profiles  for  strings  SI  and  S2  are 
quite  different,  as  shown  in  Fig.  4.  These  obstacles  profiles  have 
been  traced  from  the  most  affected  points  of  both  PV  generators: 
middle  point  of  lower  border  for  SI  and  middle  point  of  upper 
border  for  S2. 

These  two  strings  have  been  chosen  because  of  their  similarity: 
same  number  of  PV  modules  and  connection  to  the  same  inverter 
model.  The  PV  generators'  peak  power  at  standard  test  conditions 
was  characterized  at  the  beginning  of  this  study  according  to  the 
procedure  described  by  Caamano  et  al.  [37].  Since  the  strings  are 
connected  to  the  same  inverter  model,  the  MPPT  strategy  in  both 
cases  is  the  same.  It  is  important  to  note  that  the  MPPT  efficiency 
of  the  inverters  (and  its  associated  error)  is  still  present  in  the 
results  presented  in  the  next  section.  The  conversion  efficiency  of 
both  inverters  was  also  characterized  experimentally  by  measur¬ 
ing  simultaneously  inverters'  input  and  output  during  a  day  with 
clear  sky  conditions,  therefore,  the  characteristic  parameters 


1  It  must  be  noted  that  the  international  standard  describing  the  requirements 
for  the  measurement  of  irradiance  has  changed  since  the  deployment  of  this 
system  and  that  the  irradiance  sensor  does  not  follow  the  requirements  of  the  new 
version  of  this  standard.  In  particular,  the  voltage  of  the  cell  used  as  sensor,  which  is 

polarized  by  a  shunt  resistor  that  allows  measuring  the  current,  sometimes  exceeds 
the  limit  of  3%  of  V0c  as  indicated  by  the  new  version  of  the  standard.  This  leads  to 
a  slight  underestimation  at  high  values  of  incoming  irradiance.  However,  this  effect 
is  small  in  crystalline  silicon  solar  cells.  This  is  corroborated  by  the  fact  that  the 
errors  of  this  methodology  at  high  irradiance  are  lower  than  the  errors  at  low 
irradiance  (see  Section  4). 


obtained  reflect  inverters'  real  operation  conditions  (PV  generator 
power  voltage  curve).  For  this  experiment  a  precision  wattmeter 
(model  LMG  500  by  ZIMMER  Electronic  Systems  GmbH)  was  used. 
The  most  relevant  parameters  of  PV  systems  SI  and  S2  are 
summarized  in  Table  1.  The  values  in  table  were  obtained  directly 
from  the  documentation  provided  by  the  manufacturers  and  the 
installer  except  those  marked  with  a  star  (*)  which  were  measured 
by  the  authors. 

3.2.  Data  set 

The  monitoring  system  records  every  minute  instantaneous 
values  of  the  following  parameters:  in-plane  irradiance,  cell 
temperature,  ambient  temperature,  DC  power  and  AC  power. 

As  mentioned  before,  the  reference  cell  used  to  measure 
irradiance  has  the  same  orientation  and  slope  as  the  PV  genera¬ 
tors.  As  the  radiation  model  described  in  Section  3.1  calculates 
irradiation  over  a  tilted  surface  from  horizontal  irradiation  it  is 
necessary  to  apply  the  criterion  from  Eq.  (3)  to  extract  values  of 
horizontal  irradiation  from  a  typical  meteorological  year  of  Madrid 
[7],  Different  values  for  the  time  interval  w  have  been  tested:  0, 1, 
5, 10, 15  and  30  days. 

The  power  on  the  AC  and  DC  sides  of  the  PV  system  has  been 
measured  by  the  built-in  monitoring  systems  of  the  inverters.  In 
order  to  determine  its  precision,  measured  power  data  were 
compared  with  measurements  recorded  by  a  precision  wattmeter. 
For  both  inverters  the  mean  absolute  error  and  its  associate 
standard  deviation  are  given  in  Table  2. 

The  PV  system  under  study  has  been  monitored  during  one 
year:  from  March  2008  to  February  2009.  Before  analysing  model 
performance  all  records  were  checked  in  order  to  detect  faulty 
data.  Five  days  were  excluded  from  the  analysis  because  of 
partially  or  completely  missing  data.  Additionally,  11  days  were 
also  excluded  from  the  analysis  in  PV  system  SI  because  the 
inverter  monitoring  system  was  not  operating  properly.  Total 
number  of  excluded  days  amounts  to  less  than  5%  of  the  complete 
data  set 


4.  Results  and  discussion 

The  methodology  described  in  Section  2  has  been  used  to 
estimate  the  electricity  produced  by  the  PV  systems  described  in 
Section  3.  The  instant  values  recorded  by  the  monitoring  system 
were  processed  in  order  to  obtain  average  values  every  10  min. 
From  these  averages  energy  produced  in  spaces  of  10  min  has  been 
calculated.  The  choice  of  a  10  min  interval  is  a  compromise 
between  precision  and  interval  size  [38,39], 

Ten  minute  estimates  and  measured  values  were  in  turn 
averaged  to  get  hourly  values.  A  further  filtering  process  was 
applied  in  order  to  remove  values  corresponding  to  in-plane 
irradiation  lower  than  5  Wh/m2.  This  filtering  avoids  the  high 
uncertainty  related  to  low  values  of  power  but  it  was  kept  as  low 
as  possible  so  that  validation  of  the  proposed  methodology 
includes  all  weather  conditions. 

Comparison  between  the  results  of  the  proposed  methodology 
and  real  values  has  been  carried  out  at  hourly,  daily  and  yearly 
scales.  The  accuracy  of  daily  and  hourly  electricity  production 
estimates  has  been  determined  with  the  symmetric  mean  absolute 
percent  error  SMAPE  [40]: 

SMAPE(%)  =  (\/n)  ■  \00  ■  ZiWXi-yiW/iXi+yd  (17) 

where  n  is  the  total  number  of  points  in  the  data  set,  days  or 
hours;  x,-  is  the  measured  value  and  y,  is  the  estimated  value. 

The  SMAPE  parameter  overcomes  the  problem  found  at  early 
morning  or  late  afternoon,  when  irradiation  values  are  low  and  a 
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Fig.  4.  Obstacle  profiles  for  string  SI,  in  black,  and  string  S2,  in  light  grey. 
Horizontal  axis  is  azimuth  and  vertical  axis  is  elevation. 


Table  1 

Relevant  parameters  for  systems  SI  and  S2. 


Parameter  name 

System  1 

System  2 

Generator 

Slope 

34° 

34 

Orientation 

-8 

-8° 

Peak  power* 

2376  W 

2263  W 

Modules  in  series 

26 

26 

Modules  in  parallel 

1 

1 

Module 

NOCT 

47  °C 

47  °C 

p 

-2.3  mV/°C 

—  2.3  mV/"C 

Y 

0.004  C_1 

0.004  C-1 

Cells  in  series 

33 

33 

Cells  in  parallel 

2 

2 

Inverter 

P max, I 

2500  W 

2500  W 

kO* 

1.6% 

1.2% 

kl* 

3.0% 

2.8% 

k2* 

1.4% 

2.3% 

Table  2 

Mean  absolute  errors 

and  standard  deviation  of  inverter  built-in 

monitoring 

system. 

MAPE  (%) 

a{%) 

Inverter  SI 

1.5 

3.4 

Inverter  S2 

2.3 

3.5 

small  error  of  few  watts  in  electricity  estimation  is  overweighted 
by  a  more  traditional  accuracy  estimator  like  the  mean  absolute 
percentage  error  (MAPE).  Another  advantage  of  SMAPE  is  that  it  is 
less  sensitive  than  MAPE  to  errors  in  measured  values  (xf).  For 
annual  estimates  of  electricity  production,  since  they  are  single 
values,  the  percent  error  <5  has  been  used: 

<5=100-(yi-x,)/xi  (18) 

Additionally,  an  external  tool  has  been  included  in  the  study 
as  a  reference  for  a  better  assessment  of  the  methodology 
performance:  the  commercial  PVSyst  software  has  been  used. 
PVSyst  has  been  chosen  because  it  is  a  state  of  the  art  software 
widely  used  to  analyze  PV  systems  performance.  It  allows  a 
thorough  analysis  of  near  shadings  and  PV  strings  can  be  split  so 
that  the  behaviour  of  every  module  when  shaded  is  modelled 
individually. 

Tables  3  and  4  summarize  the  performance  in  the  estimation  of 
AC  power  by  the  models  described  in  Section  2  and  for  all  values  of 
w  used  to  apply  Eq.  (3)  as  well  as  results  provided  by  the  reference 


Table  3 

Performance  in  the  estimation  of  AC  electricity  for  system  SI  for  different  models 
and  different  window  sizes. 


Model 

(Eac) 

window  size 
(days) 

^(^AC,a) 

(%) 

SMAPEd 

(%) 

ad 

(%) 

SMAPEh 

(%) 

ah 

(%) 

Osterwald 

0 

-6.1 

8.4 

9.7 

23.6 

33.4 

Osterwald 

1 

-4.1 

6.4 

8.1 

21.4 

32.5 

Osterwald 

5 

-3.9 

6.5 

8.5 

20.8 

32.7 

Osterwald 

10 

-4.3 

7.0 

8.6 

21.1 

33.1 

Osterwald 

15 

-4.3 

6.9 

7.9 

20.8 

32.6 

Osterwald 

30 

-4.4 

7.3 

9.9 

21.3 

32.9 

Beyer 

0 

-0.9 

6.9 

7.9 

23.5 

35.3 

Beyer 

1 

-0.8 

5.1 

6.6 

20.8 

33.5 

Beyer 

5 

-0.5 

3.3 

4.1 

19.1 

32.8 

Beyer 

10 

-0.6 

3.0 

3.4 

18.9 

32.7 

Beyer 

15 

-0.6 

2.8 

2.8 

18.8 

32.6 

Beyer 

30 

-0.6 

3.1 

3.7 

18.8 

32.3 

Reference 

- 

-11.3 

11.7 

2.7 

54.7 

61.7 

Table  4 

Performance  in  the  estimation  of  AC  electricity  for  system  S2  for  different  models 
and  different  window  sizes. 


Model 

(£ac) 

window  size 
(days) 

^(^AC,a) 

(%) 

SMAPEd 

(%) 

ad 

(%) 

SMAPEh 

(%) 

ah 

(%) 

Osterwald 

0 

-6.2 

15.4 

18.7 

48.4 

52.7 

Osterwald 

1 

-3.1 

11.2 

14.7 

44.0 

51.6 

Osterwald 

5 

-2.5 

7.6 

9.0 

40.8 

50.9 

Osterwald 

10 

-3.3 

6.7 

7.8 

40.6 

51.2 

Osterwald 

15 

-3.3 

6.3 

6.7 

39.9 

51.1 

Osterwald 

30 

-3.6 

6.6 

9.3 

39.6 

51.2 

Beyer 

0 

1.0 

14.5 

17.8 

50.1 

54.5 

Beyer 

1 

1.6 

10.5 

15.7 

47.4 

55.8 

Beyer 

5 

1.7 

7.1 

11.1 

44.4 

54.5 

Beyer 

10 

1.5 

6.4 

10.0 

44.5 

54.8 

Beyer 

15 

1.6 

5.9 

8.7 

44.3 

55.2 

Beyer 

30 

1.5 

5.7 

7.2 

43.7 

54.8 

Reference 

- 

-14.1 

15.1 

5.9 

69.4 

67.8 

Kt 

Fig.  5.  Error  in  daily  electricity  estimations  versus  clearness  index,  comparison 
between  Osterwald  model  and  Beyer  model  for  system  SI. 


model  (PVSyst  tool).  Table  3  shows  the  results  for  system  SI  and 
Table  4  for  system  S2.  For  daily  and  hourly  values,  along  with  the 
average  error  (SMAPE)  its  associated  standard  deviation,  a,  is  also 
given.  Production  estimates  for  system  SI  are  more  accurate  than 


D.  Masa-Bote,  E.  Caamafio-Martin  /  Renewable  and  Sustainable  Energy  Reviews  31  (2014)  492-500 


499 


Table  5 

Average  error  (SMAPE)  in  daily  electricity  estimation  for  system  SI. 


Model 

All  (%) 

Sunny  (%) 

Partially  overcast  (%) 

Fully  overcast  (%) 

Osterwald 

6.5 

2.8 

5.4 

19.0 

Beyer 

3.2 

1.9 

2.6 

8.8 

Reference 

11.7 

12.3 

11.6 

10.2 

Table  6 

Average  error  (MAPE)  in  daily  electricity  estimation  for  system  S2. 

Model 

All  (%) 

Sunny  (%) 

Partially  overcast  (%) 

Fully  overcast  (%) 

Osterwald 

9.0 

5.1 

7.5 

17.6 

Beyer 

7.1 

4.5 

7.6 

16.6 

Reference 

15.1 

15.8 

15.2 

13.5 

Table  7 

Annual  irradiation  and  daily  mean  irradiation  for  every  day  type. 


Irradiation 

All 

Sunny 

Partially  overcast 

Fully  overcast 

(kWh/m2) 

(kWh/m2) 

(kWh/m2) 

(kWh/m2) 

Total 

1641.5 

814.2 

732.6 

82.4 

Daily  (mean) 

4.5 

6.1 

4.5 

1.3 

those  for  system  S2  in  all  cases.  This  result  is  not  unexpected 
because  system  S2  has  a  more  severe  obstacle  profile  than  SI. 

Comparison  of  results  for  different  window  sizes  shows  that  a 
value  of  0  for  w,  performs  significantly  worse  than  any  other  value 
for  both  systems;  the  remaining  w  values  exhibit  a  similar 
performance  with  absolute  differences  smaller  than  2%  within 
models.  For  annual  estimations  a  value  of  5  days  for  w  provides  the 
best  results  in  both  systems.  For  daily  and  hourly  values  the  lowest 
error  corresponds  to  w  equal  to  5,  10  or  15  days,  depending  on 
model  and  system.  The  differences  in  performance  for  these  three 
values  are  very  close,  with  an  absolute  difference  smaller  than 
0.5%.  However,  the  lowest  standard  deviation  error  is  achieved 
with  w=15  for  all  cases.  These  results  show  that  windows 
corresponding  to  values  of  0  and  1  for  w  are  small  in  size  and 
provide  poorly  matching  days.  A  bigger  value  for  w  (5,  10  or  15) 
performs  better  since,  for  each  day,  a  wider  collection  of  days  exist 
from  which  to  select  the  most  adequate  day.  Further  increasing  the 
value  of  w  does  not  increase  model  performance  anymore;  a  value 
of  w  =30  does  not  outperform  smaller  values,  it  only  improves 
performance  slightly  for  system  S2. 

Regarding  comparison  between  models,  the  Beyer  model 
provides  the  best  results  for  system  SI  while  for  system  S2  Beyer 
model  performs  better  for  annual  and  daily  values  while  Oster- 
wald  model  performs  better  for  hourly  values.  However,  the 
differences  between  these  two  models  are  very  narrow  for  both 
systems.  When  comparing  with  the  reference  model,  both  models 
tested  in  the  study  perform  better.  Interestingly,  the  reference 
model  performs  comparatively  worse  in  the  estimation  of  annual 
energy,  the  annual  error  being  an  order  of  magnitude  larger.  In 
addition,  this  model  underestimates  electricity  production  as 
annual  error  is  negative  for  both  systems.  This  fact  is  corroborated 
by  the  mean  biased  error  (MBE),  not  shown  on  the  tables,  which  is 
—  11.1%  for  system  SI  and  -13.6%  for  S2. 

Fig.  5  shows  the  errors  in  daily  electricity  estimation  versus 
clearness  index  for  system  SI  and  Beyer  and  Osterwald  models. 
Daily  errors  are  bigger  for  small  values  of  clearness  index  for 
Osterwald  model  while  for  Beyer  model  daily  error  is  uncorrelated 
with  respect  to  clearness  index.  These  results  corroborate  the 


appropriateness  to  include  a  model  that  incorporates  the  effect  of 
irradiance  level  on  PV  module  efficiency. 

Finally,  the  performance  of  the  proposed  methodology  has 
been  compared  with  regard  to  meteorological  conditions.  The 
365  days  analyzed  were  divided  into  three  categories: 

-  Clear  days:  the  sky  is  free  from  clouds,  diffuse  fraction  is  below 
10%  and  the  irradiance  profile  for  the  whole  day  has  a  very 
characteristic  bell-like  shape. 

-  Partially  casts  days:  the  presence  of  scattered  clouds  in  the  sky 
covers  and  uncovers  the  sun  repeatedly  and  irradiance  changes 
its  value  very  fast. 

-  Fully  overcast  days:  the  sky  is  totally  covered  by  clouds; 
irradiance  is  irregular  over  the  day  but  with  low  values,  never 
exceeding  600  W/m2. 

Of  the  365  analyzed  days  analyzed  136  were  classified  as  sunny, 
164  as  partially  overcast  and  62  as  fully  overcast.  Three  days  could 
not  be  classified  because  of  missing  data. 

Average  daily  electricity  estimation  errors  (SMAPE),  separated 
into  day  types  are  shown  in  Tables  5  and  6  for  Osterwald  and 
Beyer  models,  using  w=5  days  as  half  window  size.  Average  daily 
errors  obtained  from  the  reference  model  are  also  shown.  The 
previous  trends  observed  for  all  days  also  apply:  electricity 
estimation  is  more  accurate  for  system  SI  than  for  S2  and  the 
Beyer  model  performs  better  than  Osterwald  model.  For  these  two 
models,  performance  for  sunny  and  partially  overcast  days  is 
significantly  better  than  for  fully  overcast  days,  being  performance 
on  sunny  days  the  best  overall.  The  reference  model,  on  the  other 
hand,  exhibits  the  opposite  behaviour:  it  behaves  better  on  fully 
overcast  days  than  on  sunny  and  partially  overcast  days.  For  fully 
overcast  days,  the  reference  model  performs  better  than  any  other 
model  for  PV  system  S2  and  it  is  comparable  (slightly  worse)  to 
Beyer  model  for  PV  system  SI.  It  must  be  noted  that  the  impact  of 
fully  cast  days  in  annual  electricity  production  is  significantly  less 
than  for  clear  and  partially  overcast  days  because  of  the  differ¬ 
ences  in  irradiation  for  every  day  type.  In  this  sense,  annual  global 
irradiation  and  its  breaking  down  into  day  typologies  are  pre¬ 
sented  in  Table  7.  Mean  values  of  daily  irradiation  have  been  also 
included  to  enable  comparisons  between  different  weather  con¬ 
ditions  independent  of  the  number  of  occurrences  of  every  day 
during  the  year.  Average  daily  irradiation  for  partially  overcast 
days  is  74%  of  irradiation  available  on  sunny  days  while  the 
percentage  for  fully  overcast  days  is  as  low  as  21%.  When  looking 
at  accumulated  irradiation  during  the  year,  50%  of  total  annual 
irradiation  corresponds  to  sunny  days,  45%  to  partially  overcast 
days  and  only  5%  to  fully  overcast  days  while  17%  of  days,  62  out  of 
365,  were  classified  as  fully  overcast. 

5.  Conclusions 

A  new  methodology  for  estimating  the  electricity  produced  by 
BIPV  systems  has  been  presented.  This  methodology  incorporates 
a  novel  model  for  shading  losses  estimation  developed  specifically 
for  BIPV  systems.  This  model  for  shading  losses  estimation  is  based 
solely  on  geometrical  considerations,  an  electrical  model  for  the 
PV  generator  or  information  about  modules  internal  structure  and 
wiring  are  not  needed.  The  methodology  is  validated  with  models 
to  estimate  the  useful  electricity  a  PV  system  can  produce.  Two 
alternative  models  to  estimate  electricity  production  have  been 
tested.  Noteworthy,  the  model  that  estimates  electricity  without  a 
direct  measurement  of  PV  generator  peak  power,  not  usually 
performed  on  BIPV  systems,  provides  results  comparable,  better 
in  some  cases,  than  the  model  requiring  an  explicit  measurement 
of  such  characteristic. 
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The  methodology  has  been  evaluated  during  one  year  with  real 
data  from  two  PV  systems  representative  of  BIPV.  The  evaluation 
has  been  performed  over  a  wide  range  of  irradiance  values, 
starting  from  values  as  low  as  5  W/m1 2 3 4 5 6 7 8 9 10 11 12,  covering  all  weather 
conditions.  The  error  obtained  in  annual  electricity  estimation  is 
below  1%  for  the  best  performing  model.  The  SMAPE  of  daily 
electricity  is  below  3%  for  the  same  model.  It  has  been  also  found 
that  including  the  effect  of  reduced  performance  at  low  irradiance 
levels  improves  the  estimation  of  the  electricity  yielded  by  a  PV 
system.  When  compared  to  a  reference  model  widely  used  in  PV 
system  performance  analysis  that  can  be  considered  representa¬ 
tive  of  the  actual  state  of  the  art,  the  methodology  presented 
performed  better  than  in  annual,  daily  and  hourly  estimation  of 
electricity  production. 

The  proposed  methodology  is  simple,  easy-to-use  and  can 
provide  fast,  accurate  results  at  low  costs  since  it  does  not  require 
specific  equipment.  These  features  make  this  technology  suitable 
for  application  in  the  supervision  of  BIPV  systems.  Although  the 
methodology  has  been  designed  specifically  for  building  inte¬ 
grated  PV  systems  it  can  be  applied  also  to  other  types  of  PV 
systems  different  than  building  integrated. 
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